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The gyrocenter phase-space transformation used to describe nonhnear gyrokinetic theory is redis- 
covered by a recursive solution of the Hamihonian dynamics associated with the perturbed guiding- 
center Vlasov operator. The present work clarifies the relation between the derivation of the gyro- 

' center phase-space coordinates by the guiding-center recursive Vlasov method and the method of 

' Lie-transform phase-space transformations. 
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I— i' is to expand the Vlasov operator L = d/dt = Lq -\- e Li + L2 + ■ ■ ■ and the Vlasov distribution / = /o + e /i + • • • 
i-C , asymptotically in powers of a small ordering parameter e. One then proceeds with a recursive solution of the Vlasov 
^p-i' equation ([T]) at each order in the hierarchy: 

S : = Lo/o 1 

^ . = Li/o + io/i 

riS ■ = L2/0 + ii/i + Lof2 } ■ (2) 

Oh! 

Two important assumptions are associated with the hierarchy ([2]). First, we assume that the lowest-order equation 
^ Lofo = has a known (exact) solution. For example, the lowest-order dynamics characterized by Lq is often associated 
with a cyclic (lowest-order) orbital angle ip (i.e., Lq = uj^d/dip, where uj^ = dip/dt) so that the lowest-order Vlasov 
equation Lq/o = simply implies that /o is independent of the orbital angle ip. Second, we assume that the operator 
Lq can be inverted (e.g., L^^f = u)~^ J f dip) so that the solution for the first-order correction /i in ([2]) may be written 
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as /i = /i — L^^(Li/o), where fi satisfies the homogeneous equation Lo/i = 0. The small ordering parameter e 
appearing in this asymptotic expansion is defined by the relation L^^ Lk = 0{e^). 
I The purpose of the present paper is to compare and contrast the derivation of gyrocenter phase-space coordinates 
\ by the guiding-center recursive Vlasov method and the method of Lie-transform phase-space transformations. In §|TT] 

■ the expansion of the Vlasov operator L — d/dt defined in ([1]) is given in powers of a small ordering parameter e. The 
\ guiding-center recursive Vlasov (gcrV) method is introduced in ij lllll based on expansion of the guiding-center Vlasov 
. operator Lgc = T^^^LTgc defined in terms of the guiding-center push-forward (T~(.^) and pull-back (Tgc) operators. 

■ In order to keep the analysis focussed on gyrokinetic applications, we assume that the background magnetic field 
0^ , is uniform and that the particles move under the influence of a fluctuating electrostatic field (with slow and fast 

' space-time scales). Using the gcrV method, we derive explicit expressions for the gyrocenter phase-space coordinates. 
I> ! In ^ llVi the gyrocenter Lie- Transform (gyLt) method is applied to the derivation of the gyrocenter phase-space 
coordinates. In order to allow comparison with expressions derived by the gcrV method, higher-order terms are kept. 
In §|Vl the gyrokinetic Vlasov equation is derived both by the gcrV and gyLt methods. While the gcrV derivation 



^ . does not allow for a systematic truncation scheme that preserves energy conservation, the gyLt method does since 
■ ■ ' it is naturally associated with a variational formulation. We also introduce the gyrocenter pull-back operator Tgy 
and discuss its physical interpretation. In i) IVI[ we present the gyrokinetic Poisson equation expressed in terms of the 
gyrocenter moment (with respect to the gyrocenter Vlasov distribution) of the gyrocenter push- forward 'T~y5'^^ of the 
guiding-center delta function ^g^. = Tgp^5^(x — r). Lastly, our work is summarized in ij lVIII and Appendix A presents 
the guiding-center phase-space transformation for nonuniform magnetic fields. 

II. EXPANSION OF THE VLASOV OPERATOR 



We begin with the expansion of the Vlasov operator 

1 = 1 + (.„b + v,) .V + ^(e + Xxb) 4 . Lo + eLr, (3) 
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which is written in terms of the two operators Lq and Li (defined below). Here, within the context of the derivation of 
gyrocenter coordinates for an electrostatic gyrokinetic model , we consider the case of a strong background magnetic 
field and a small quasi-static electric field E = — Wcj). To simplify the analysis, we also consider a uniform background 
magnetic field (i.e, the unperturbed guiding-center dynamics is represented simply in terms of parallel motion along 
straight magnetic field lines). 



A. Zeroth-order Vlasov Operator 



In a uniform magnetic field, the zeroth-order Vlasov operator Lq is defined as 

d 



(4) 



where G is an arbitrary function on particle phase space and the perpendicular velocity vj^ = fldpg^/dO is expressed 

in terms of the gyroradius vector p^^ — (b/fi) X v, which depends explicitly on the gyroangle 6. In the second 
expression in (jj]), we introduced the guiding-center pull-back Tgc = exp(— Pg^, • V) and the guiding-center push- 
forward = exp(pgj.'V), which are both associated with the guiding-center phase-space transformation for a 
uniform magnetic field (the guiding-center pull-back and push-forward operators are given for a nonuniform magnetic 
field in Appendix A). Note that the second term (vj^ • V) in ([4]) may be of the same order of magnitude as the first 
term [ild/dO) when it is applied to short-wavelength fluctuating fields (e.g., fluctuations that satisfy the gyrokinetic 
ordering j|]). 

After defining the guiding-center push-forward of an arbitrary function G on particle phase space: 

Ggc = Tg^^G, (5) 
we arrive at the final expression for the zeroth-order Vlasov operator: 

LqG = Tgc qq"^ ^ '^sc (-^ogcGgc ^ , (6) 
where Lggc = ^d/d9 is the lowest-order guiding-center Vlasov operator. This operator can easily be inverted: 

L-giF = n-' jFde = n-' j (f - {f)) de, (7) 

where F denotes the gyroangle-dependent part of F and (F) denotes the gyroangle-averaged part. We immediately 
see that the zeroth-order operator explicitly involves the lowest-order guiding-center transformation (for a uniform 
magnetic field) and we note that this formulation can also be applied to the bounce-motion Vlasov description of 
magnetically-trapped particles 0]. 



B. First-Order Vlasov Operator 

The first-order Vlasov operator 

is explicitly decomposed in terms of the parallel and perpendicular components of the electric field. The ordering 
Lq + e Li in ([3]) implies that the gyromotion time scale is the shortest time scale in our plasma physics problem. 
In particular, the ordering fl~^Li± ^ 1 implies that the perpendicular E x B velocity is small compared to the 
characteristic (i.e., thermal) velocity of a particle (this is the drift ordering [l3|)- The ordering n~^Lm <C 1, on 
the other hand, implies that the time scale of interest is long compared to the gyration period and that the parallel 
gradient length scale is long compared to the gyroradius. These two orderings are consistent with the guiding-center 
and gyrocenter orderings fsj. [Note that it is also possible to work with an alternate ordering where ri~"'^Li|| ^ n~^Li^, 
in which case one could elevate Li|| = L2 to a higher order in ^.] 
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We first look at the operator Li± in ([8]), which can be rewritten as 

'9pg, dG dpg, dG 



Li±G 



g- El • 



B 



89 dn dfi 36 



B 
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dp^A dG 



89 J 8/1 



V4> 



dpgc\ dG 



8fi J 89 



(9) 



where /i = m\'v±\'^ /2B denotes the guiding-center magnetic moment (a guiding-center invariant in a uniform magnetic 
field) and p denotes the gyroradius vector. We now use the identity 



dp 



8J' 



T, 



8J' 



where 9/9 J* = {8/89, 8/8y) and i/igc 
operator ([9]) can be written as 



T; 



(X + p ) defines the guiding- center scalar potential, so that the 
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89 J 8^i 



^ - T, 



90, 
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8^1 J 89 



gc 



TgcGgc, 



where G = TgcGgc is expressed in terms of the pull-back of Ggc- Next, we note that the operators 8/8,P and 
Tgc = exp(— • V) do not commute: 
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^TgcGgc ^ 



/8G. 



\ 8 J' 



Hence, using the fact that the pull-back operator is distributive Tgc(i^ G) = (TgcF) (TgcG), we obtain 



Li±G 



n 

B 



9- Tgc 



"gc 



8G, 



8G, 



gc 



89 8ii 8ii 89 

fdPgc 9pgc 9pgc 9pg 
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\ 89 9/1 9/i 89 



•VGh 



(10) 



where we have used the explicit property 90gc/9J' — V(/)gc ■ 8p^^/8,P to obtain the last two terms. Lastly, by using 
the identity (valid for any vectors F and G) 



f8p 8p„^ 8p„^ 8p„ 



V 89 8^l 8^1 89 
we obtain our final expression for Li±^: 



•G 
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(FxG), 



ii±G 



9T 



bgc 9Ggc 



9(/)gc 9Gg 



B \ 89 8^1 9/1 961 

9 Tgc gc, Ggc}^ 
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where we introduced the "perpendicular" components of the guiding-center Poisson bracket. 
We can similarly write ii||G in terms of 0gc and Ggc as 



I/i||G = Tgc ( ^gc + 9 {Ggc, 



gc / ' 



(11) 



(12) 



where dgc/dt = 8/8t + ti|| b • V is the guiding-center Vlasov operator (in a uniform magnetic field), and 



dF 



l|gc 



dvn 



8vii 



denotes the "parallel" guiding-center Poisson bracket. By combining the perpendicular and parallel components (jlip 
and l|12p. we therefore obtain the final expression for the first-order Vlasov operator 



LiG - Tgc ( — Ggc + q {Ggc, 4>gc}gc ] = Tgc ( iigcG 



dt 



(13) 
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where { , }gc now denotes the full guiding-center Poisson bracket. We note that the guiding-center representation 
(|13p also holds for a nonuniform background magnetic field, where the guiding-center Poisson bracket now contains 
corrections associated with magnetic-field nonuniformity (see Appendix A) . 

Lastly we note that it was the original insight of Catto [4] to recognize that the transformation from particle 
coordinates to guiding-center coordinates greatly simplifies the recursive solution of the Vlasov equation ([T]) within 
the context of linear gyrokinetic theory. What we have shown in this Section is that this simplification naturally 
extends to the inclusion of the guiding-center Poisson bracket { , }gc as well as finite-Larmor-radius (FLR) effects 
associated with the scalar potential (j) (through the pull-back and push-forward operators) in the first-order guiding- 
center Vlasov operator -/^igc. 



III. GUIDING-CENTER RECURSIVE VLASOV METHOD 



In this Section, we use the guiding-center recursive Vlasov (gcrV) method to derive asymptotic expansions for the 
gyrocenter phase-space coordinates. Here, the gcrV method is defined in terms of recursive solutions involving the 
guiding-center Vlasov operator 



T = T^^ T T 
^gc — ' gc ' g 



de 



dt 



"gcjgc 



gc 



L 



Igc- 



(14) 



The gcrV method therefore naturally takes into account the full FLR effects of the perturbation scalar potential. 
A generic gyrocenter variable Qgy can be expressed in terms of a time-dependent asymptotic expansion 



Qg 



VgV 



gy 



— Tgy Qgc — Tgy (TgcQ) 



(15) 
Vgy 

is decomposed in terms of gyroangle-independent and gyroangle-dependent parts (respectively). The asymptotic 
expansion (jlSp is constructed by gcrV method by requiring that the quantity Qgy = igcQgy is gyrophase independent. 
This condition yields the following fcth-order expressions 



where Qgy'' = Qgc denotes the guiding-center variable and the fcth-order gyrocenter variable Qgy' = 



'gy 



'gy 



-^lgc^;gy 



(Llgc) Qgy 



iigcQ^'^ 



and 



Wgy — -^Ogc 



-'^IgC^fgy 



(16) 



(17) 



Note that the solution for the gyroangle-independent part Q^^ appears at the fcth order. 

The gcrV method now proceeds as follows. At zeroth-order, for the gyrocenter phase-space variable Zq^^ ^ 9, we 
find the definition 



i^gy 



= n 



Ogy 



de 



= 0, 



(18) 



so that the zeroth-order gyrocenter (guiding-center) coordinate Z^^^ ^ 9 must be independent of the gyroangle 6 (i.e., 

(Z^gy) = Z^gy). Obviously, (^gy)o = ^ for the gyrocenter gyroangle Z^^y 0. 
At first order, we find the gyroangle-independent expression 



(^gy)^ = (^^Sc) ^Ogy 



"gc^Ogy 

dt 



1 {^Ogy 



o ^ "gy^Ogy 
''>SC - dt ' 



(19) 



where dgy/dt denotes the lowest-order gyrocenter Vlasov operator (which includes the E x B velocity and parallel 
electric field associated with ((/)gc)). The gyroangle-dependent first-order expression, on the other hand, is obtained 
from P?|) as 



-^Ig 



L, 



Ogc 



l^Ogyi '?^gc| 



\ ^Ogv *^gc f ; 
I. J gc 



where 



i^gc 



d9 = 



i^gc 



) d9 



(20) 



(21) 



5 



denotes the indefinite gyroangle integral of the gyroangle-dependent part of (/)gc. Note that the gyroangle-independent 

part of Zfgy (denoted ^igy = ^fgy ~ ^igy) must be determined at the second order. 
At second order, we find the gyroangle-independent expression 

and the gyroangle-dependent expression 

Z2gy = ~ ^ {^Igy ^Sc}^^ - -^Ogc (-^Igc -^fgy) • (23) 

In ([22|) . we see that the role of ^igy is to ensure that (^gy)2 satisfies any desired property we want (e.g., be gyroangle- 
independent or zero). One clearly sees how the gcrV method can be extended to higher order. 



A. Gyrocenter magnetic moment 

The easiest gyrocenter phase-space variable one can construct by the gcrV method is the gyrocenter magnetic 
moment figy since it is to be constructed as an invariant of the gyrocenter dynamics, i.e., /igy = to arbitrary order 
in e. At zeroth order, we easily find iogcM = Oj where /^ogy = M denotes the guiding-center magnetic moment (which 
satisfies dgc/^/rfi = 0). 

At first order, the requirement (/igy)i = becomes 

which is easily solved as 

Mlgy = Mlgy + ^^gc ^ {THsy ~ ^ ('/'gc)) + -| 0gc, (25) 

where /J^gy denotes the gyroangle-independent part of the first-order gyrocenter magnetic moment (to be determined 
at the second order). We note that the first-order magnetic- moment correction Afj, derived by Taylor [l^ can be 
expressed as the guiding-center pull-back A/i = Tgc Jj-igy of the gyrocenter magnetic-moment correction (j25p . where 

Tgc^gc = Tgc(Tg^V - {4>gc)) ^ (j) - Tgc((?igc). (26) 
At second order, the requirement (/igy)2 = becomes 

= + {Mlgy, '^gclg. + ^ (27) 

rfgc /_ , ~ ^ , „ / /^T7 lA \\ A, \ in ^/^2gy 



(Mlgy + A^lgy) + 9 {(^Igy - ;|('/'gc)) , (/-gc} ^ + ^ 



l/^lgy ■ ^^^y, • ^ ly^^gy ^ ^^g., j , ^g. j ■ - 

where we used (|25|) for ^igy (with {4>gc, 0gc}gc = 0) and the role of /Iigy is to ensure that the gyroangle-independent 
right side of (|27p is zero. This condition yields the homogeneous equation dgyjli^y / dt = 0, whose solution is simply 
/J^gy = 0. The solution to the gyroangle-dependent part of ([27| is finally expressed as 

_ q dgyigc . . 

M2gy = M2gy " ^ , (28) 



where 



dt ydt • B"-^-"^^' -y^g-^-^^s i^sc 

and the gyroangle-independent part 7l2gy is determined at the third order. 
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Lastly, it is instructive to compute the third-order component of the gyrocenter magnetic moment. Here, (/igy)3 = 
becomes 

_ i^gcM2gy , „ r„ ,1 , o ^M3gy 

= ^ (M2gy + M2gy) + q { (Msgy + M2gy) , 0gc}^^ + ^ . (29) 

The role of 'p2gy is to ensure that the gyroangle- independent right side of (|29p is zero, which implies that 7l2gy is a 
solution of the inhomogeneous equation 



'^gyM2gy rT_\\_ 1^ / } 2 '^gy^gc I 

Jgc/ 



Because the right side of this equation is non- vanishing, there must be a non-trivial solution for 7l2gy By using 
properties of the guiding-center Poisson bracket { , }gc and the gyrocenter Vlasov operator dgy/dt (after tedious 
algebra), we find the second-order solution 

This solution is more trivially obtained with the Lie-transform approach presented in ij llVl To lowest FLR order, 
where (figc = d^gc/d6 ~ • Vj_0, we note that 



A^2gy ^ ({*gc. ^,c}^)- (30) 



{$gc, ^sc}^ 

so that 



» g((0gc)^) ^ ^ \y±4>\' 

B 9/i m n 



M2gy - ^ |UEP, (31) 



where ue = (cb/B) X denotes the E x B velocity. We omit the explicit derivation of Jlsgy, which is obtained 
from the gyroangle-dependent of the third-order equation (|29p and contains terms that are of second order in 
Up to second order in e (and first order in 51^^), the gyrocenter magnetic moment is therefore expressed as 




.„ = . + %l*.c-^4|^ ). (32) 



We note that, to lowest FLR order, the first two terms in (f32|) appear naturally in the expansion of the magnetic 
moment figy = m |vj_ — (u^ -|- Up)p/2i3, where up = — {c/ BQ) dVcf)/ dt denotes the polarization drift velocity. 
Lastly, we note that Parra & Catto [ll| only computed the first-order correction {q(j>gc/B) to the gyrocenter magnetic 
moment. While the polarization-drift correction dgy^gc /dt is generally not kept in standard gyrokinetic theory 
we show in ij llVI how it also appears naturally in the Lie-transform approach. 



B. Gyrocenter gyroangle 

At zeroth order, we easily find Logc^ = where ^ogy = ^ denotes the guiding-center gyroangle. At first order, we 
find 



which yields the gyroangle-independent equation 

while the gyroangle-dependent equation yields the solution 

^Igy = - = 6lgy - Olgy, (35) 

where 9igy denotes the gyroangle-independent part to be determined at second order. 
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C. Gyrocenter position 

The zeroth-order gyrocenter position is the guiding-center position Xogy = X = r — p^^, where r is the particle's 
position and p^^ is the gyroangle-dependent gyroradius vector. At zeroth-order, we easily find LogcX = 0, since the 
guiding-center position X is independent oi 6. At first order, we find 



rfprX , , ^ „ 9X 



igy 



'gc QQ 

= t;||b + ^xV^gc + (36) 
which yields the gyroangle-independent first-order equation for the gyrocenter velocity 



Xgy)^ = «||b + ^xV(0gc). (37) 

We therefore see that the lowest-order gyrocenter motion is described in terms of parallel motion along and E x B 
motion across the field lines. The gyroangle-dependent equation obtained from p6p yields the solution 

~ cb ~ — 

^Igy = ^ X ^*gC = Xlgy - Xlgy, (38) 

where Xigy denotes the gyroangle-independent part to be determined at second order. 

D. Gyrocenter parallel momentum 

The zeroth-order gyrocenter parallel momentum is the guiding-center parallel momentum P||ogy = At first 

order, we find 

(P||gy)l = -9b.V0ge + (39) 
The gyroangle-independent part of this equation yields 

(P||gy)i - -gb-V(</>ge), (40) 

while the gyroangle-dependent part yields 

1^ 

Plligy = ^^-V^-gc = Piiigy " P||igy, (41) 
where P\\igy denotes the gyroangle-independent part to be determined at second order. 

E. Gyrocenter kinetic energy 

We now use the operators Logc and iigc to derive an asymptotic expansion for the kinetic energy Kgy = Kogy + 
e Kigy + ■ ■ ■, where 

Kogy = mv'^/2 = p^/2to + = Kg^ 

is expressed in terms of the lowest-order guiding-center coordinates p\\ and fi. We want to construct Kgy such that 
Kgy is gyroangle independent. At zeroth order, we easily find that {Kgy)o = LogcKogy = 0, i.e., the guiding-center 
kinetic energy Kgc is a constant on the gyromotion time scale. 
At first order, we find 

(^gy)l = + 1 {^Ogy, <^gc}g. + 

= -||b-V0g.) + (42) 
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where we used the fact that dgcKogy / dt = 0. The gyroangle-independent and gyroangle-dependent parts of (|42)) yield 



{Kgy)i = - qvn b- V( 



and 



K 



igy 



n 



Vn b-V$„ 



K 



Igy 



(43) 
(44) 



where the gyroangle-independent part i^igy is determined at the next order. 
At second order, wc find 



(^gy)2 



dgcKl^ rj^ , , O 9K2gy 



dt 



"'gy 



if 



Igy 



gc ■ -- gg 

Klgy + Klgy I , (Pgc, 



(45) 



{(' 



y) ' '/'gel 



2gy 
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By choosing if igy = in (l45l) . we obtain the gyroangle-independent part 

(^gy)2 = g (^[Kigy, ?gc}g^y 
Next, by using the Jacobi property of the guiding-center Poisson bracket, we introduce the identity 

{^gc, {$gc, G}gc}^^ = ^ {g, {$gc, 0gc}gc}^^ 

+ I {^.c, {5gc, G}ge} 



(46) 



where G is a gyroangle-independent function and (j>gc = d^gc/dO. Substituting into 
(|i7)) . the second-order gyrocenter kinetic equation becomes 



(47) 

and using the identity 



(ifgy)2 



Kgc, ^ ({^-gc, 0gc}gc 



gc 



The gyroangle-dependent part of the second-order kinetic energy, on the other hand, is expressed as 



if. 



2gy 



}_ ^gy^igy 
n dt 

q_ dgy$gc 
SI dt 



q_ 



|iflgy, 0gc}^^ - <^{^lgy, '^^gc}^^^^ 



de 



(48) 



where we have ignored terms of order fl ^ and the gyroangle-independent part if 2gy must be computed at third order. 
The gyrocenter kinetic energy ifgy is therefore expressed as 



Kgy = ifgc + g^Z-gc + b.V$gc - ^ 



Kgc + q<Pgc 1 



q d{(j)g 



B d^i 



n dt 

q I d cb 
n \dt^^ 



(49) 



X V( 



•V 



where terms of second order in ^ have been omitted and the terms of first order in 17 ^ associated with the 
gyroangle-independent part if 2gy have not been computed. While Parra & Catto [ll| captured the first-order term 

correctly, their second-order term includes — {q/n) dt^gc only and ignores the second-order correction terms due to 
((/)gc). In fact, Parra & Catto [ll[ systematically ignore (fP-teims in their derivations of gyrokinetic variables except 
in their revised quasineutrality condition. 

Lastly, the gyrocenter equation for ifgy is expressed as 



if, 



gy 



b-V(gr(0gc) - — (^{$gc, (/-gclgc 



(50) 



which includes a nonlinear (quadratic) contribution to the parallel electric field generated by {(t>gc)- 
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IV. GYROCENTER LIE-TRANSFORM APPROACH 



The transformation from the extended guiding-center coordinates Z^^ = Z" = (X,p|| , ^, 0; w, t) to the gyrocenter 
coordinates Z^y = Z ~ (X,p|| ,77, 0; w, t) is expressed as an asymptotic expansion 

T = Z- + eGt + e^i^G", + ^Cfgl) + (51) 

where the nth-order generating vector field G„ is evaluated at order e" to eliminate gyroangle-dependence in the 
Hamiltonian. Here, the energy coordinate w is introduced as the canonically-conjugate coordinate to time t. In 
Hamiltonian Lie-transform perturbation analysis (appropriate for electrostatic perturbations), the generating vector 
fields are expressed in terms of the extended guiding-center Poisson bracket { , }gc, which now includes the canonical 
pair {w, t), as 

Gl ^ {Sk, , (52) 

where the functions (5*1, S'2, ■ • •) are known as the gyrocenter gauge functions (which are assumed to be explicitly 
gyroangle-dependent). We note that the time coordinate t is unaffected by the time-dependent gyrocenter transfor- 
mation if G\. = dSk/dw = at all orders. 
The extended guiding-center Hamiltonian is 




Hgc ^ \ — + ^i,B + q{<f)gc) - w \ + e (gf(/)gc = iJogc + e i?igc, (53) 



where we have explicitly separated the gyroangle-dependent part 0gc of the guiding-center scalar potential as 
the perturbation that destroys the guiding-center mametic moment (i.e., {/x, iJogdgc = and {/i, -ffigdgc 7^ 0). 
While the separation adopted in ([55)1 is nonstandard [3|, it is consistent with the gcrV method presented in ij lHIl 
Note that this separation appears when the electrostatic potential has a large-scale component {4>gc) and a small-scale 
component (/>gc, which satisfy the ordering 

q^gc ^ T q{4>gc), (54) 

where T denotes the plasma temperature. This ordering is consistent with the generalized gyrokinetic ordering Q 
(p/Aj_) q(f> <C T, where the perpendicular gradient length scale X± scales as X± ^ p for the small-scale component (jy^c 
while it scales as A_l » p for the large-scale component ((/)gc). 

The gyrocenter transformation ([5T|) is chosen at each order so that the gyrocenter Hamiltonian 

-ffgy — Hogy + e Higy + 6^ H2gy + ' ' ' (55) 

is gyroangle-independent, where i?ogy = -^Ogc- According to Hamiltonian Lie-transform perturbation theory, the 
first-order and second-order gyrocenter Hamiltonians are 

Hlsy = Higc ~ {Si, i?Ogc}ge = -H^lgc - + ^ ~Qfj ' (^6) 



H2gy - - {Si, i?lgc}gc + 2 l*^^' ^Ogc }gc| 



where we used i?2gc = 0. It is straightforward to extend the Lie transform approach to a nonuniform magnetic field, 
since the operator dgy/dt — { ■ , i?ogc}gc is valid in general magnetic geometry, with the guiding-center Poisson 
bracket suitably generalized for nonuniform magnetic fields [see (|A8p ]. 

A. First-order analysis 

At first order, since i?igc — q4>gc and (-ffigc) = 0, the expression for the gyrocenter Hamiltonian is simply 

Higy ^ 0, (58) 
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while the first-order gauge function 5*1 is the solution of the gyroangle-dependent equation 

The reader should not be alarmed by (j58|) and remember that the gyroangle-averaged scalar potential {4>gc) appears 
in the gyroangle- independent perturbed guiding-center Hamiltonian -ffogc in (|53p . 
Up to second order in il^^, the solution for 5*1 is 



^1 - f^'fgc - ^2^^^ + ••• (60) 



where $|c^^^ = / $^c^ dO, with $|c^ = $gc and = 0gc. 

B. Second-order analysis 

At second order, the expression for the gyrocenter Hamiltonian is 

while the second-order gauge function 5*2 is the solution of the gyroangle-dependent equation 



dt 

where terms of order f2~^ were omitted on the right side of (|62|) . Up to second order in fi"^, the solution for 52 is 



^6*. (63) 



C. Gyrocenter Coordinates 

One of the main advantages of the gyrocenter Lie-transform approach is that the gyrocenter phase-space trans- 
formation is expressed solely in terms of the scalar fields (^i, iS'2, • • •)■ The gyrocenter phase-space coordinates are 
constructed from the asymptotic expansion 



Z" = Z" + e {5i, + £2 (^{^2, ^"}gc + I {Si, {Si, ^"}gc}ge) 



(64) 



2 

{S2, ^"}gc + ^ {*gc, {*gc, ^"}gc}^ 

where we substituted the expression (|5D|) for Si while the expression for 6*2 will be used only when needed. Hence, 
the gyrocenter position X is 

X = X + 6^ {$gc, X}^^ + ••• = X - eAx V$gc + (65) 
the gyrocenter parallel momentum p|| is 



{*gc, P\\} ^ + ••• = p„ e~V$gc + (66) 
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the gyrocenter gyroangle 9 is 

the gyrocenter magnetic moment Ji is 

2 _ _ 
^ ^*gc, {$gc, Migcl- 



({52, Migc + ^ {*gc, {$gc, Migc} 



= M + e ^ </>gc 



2 



^{<i>gc, <^gc}^J + (68) 



B y n dt I 2Bn 

where the expression for 5*2 was used in obtaining the last expression, and the gyrocenter energy coordinate W is 

_ _ q 5$gc 



n dt 



(69) 



The gyrocenter kinetic energy Kgy is 



n n ^ dt 

r,2 



gc 



+ (^{^2, i^gclge + ^ {*gc, {<fgc, ^gclgc}^^ 



2Vl 



B dfi J n \dt B 

2 



(70) 



where terms of second order in have been omitted. One can see that the second-order gyro-independent part 
of the kinetic energy if 2gy appears naturally in the Lie-transform approach, whereas it would require an extensive 
computation to obtain it in the guiding-center recursive Vlasov approach (we have skipped this computation, see 
i^ lIIIEl for details) . Note that the gyrocenter kinetic energy ([70)) can be exactly expressed as 

i^gy - i + J^B, (71) 

when the definitions (|66p and (|68p for p|| and /I are used. 
Lastly, the Jacobian for the gyrocenter transformation is 

Jgy = Jgc - e-g^(^Jgc {Si, ^"}gc) + ■•• ^ Jgc, (72) 
where JTgc is a constant in a uniform magnetic field. This result comes from the fact that 

^ f ^ re- 7a\ \ - ^ f -7 7^/3 _ _ -r 7a/3 ^^-^1 ^ r, 

gZ<^ i^i, ^ igc) - QZ'- [Jsc-J Q^pj - Jgc J QZ'-QzP - ^' 

which follows from the antisymmetry of the guiding-center Poisson tensor J"^ and the Liouville identities 
da{Jgc J"'^) = 0. Note that 172]) is true to all orders in e. 



V. GYROKINETIC VLASOV EQUATION 



We now derive the gyrokinetic Vlasov equation ([T]) by the gcrV method and show how the gyrocenter phase-space 
transformation (j64p generated by the Lie-transform scalar fields (5*1, S'2, • ■ •) are involved in the transformation from 
the particle Vlasov distribution / and the gyrocenter Vlasov distribution F . 
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A. Recursive Vlasov Derivation 



First, the particle Vlasov distribution / is expressed in terms of the guiding-center Vlasov distribution F by the 
push-forward operation F = T^cf — Fq + e Fi + F2 + ■ ■ where F^ = T^^f^- At zeroth, first, and second orders, 
we therefore have 



= iogc^O 

= LigcFo + LogcFi 
— LigcFi + LogcF2 



(73) 



where the guiding-center operators Logc and Ligc are defined in ([6]) and (fTS]) . 

At the zeroth order, the Vlasov equation Logc-Fb = implies that Fq is independent of the gyroangle 9. At first 
order, the gyroangle-independent part of iigc^o yields 







dt dt 

while the gyroangle-dependent part yields a solution for Fi: 



Fl ^ - ^Ogc ( Llgc Fq 



q_ 



{Fo, <fgc} 



(74) 



(75) 



The solution for the gyroangle-independent part Fi = Fi ~ Fi must come from the second-order analysis. 
At second order, the gyroangle-independent part of Lig^Fi yields 







d-gyF 1 

dt 

dgyF 1 
dt 



{Fi, 0gc} 



' 2n 



({$gc, 0gc}^J 



(76) 



while the gyroangle-dependent part yields a solution for F2: 



F2 = 



_q_ 



dt 



{Fi, 0gc} 



qi{Fi, 0gc}^^ - ({f,, 0gc} 



d0 



Fo, 



Fo, 



dt 



<1_ 



{Fu $gc} 



{{^0, $gc}^^,$g 



(^{$gc, </'gc}^^ - <^{*gc, <^gc} 



gc 



d0 



(77) 



B. Gyrocenter Pull-back Operator 



We can combine the recursive solutions ((75)) and fTTl) for the guiding-center Vlasov distribution to obtain the 
following expansion 



F = Fo + eFi + F2 



1 dgy^'gc 



= (^0 + eFi 
(2)\ ^ 



-).( 



eFi + e'F2 



e^^{^gc, F,} 



n dt 



Fn 



gc 



2n^ 



+ Fo, 



he] ^ - ({^gc, ^gc} 



de 



(78) 



gc 
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where F = Fq + e Fi + F2 + ■ ■ ■ defines the gyrocenter Vlasov distribution and we used the identity (|47| for the 
last two terms. Using the gyrocenter scalar fields {Si, S2, ■ ■ ■), the relation ([78| between the guiding-center Vlasov 
distribution F and the gyrocenter Vlasov distribution F can also be expressed as 



= F + e{Si, F}^^ + '-{Si, {Si, F}}^^ + {S2, F}^^ + ■■■ 

^ TgyF, (79) 



where the guiding-center Vlasov distribution F = Tgy F is represented as the gyrocenter pull-back of the gyrocenter 
Vlasov distribution F. We immediately see that the gcrV method has generated a solution that is exactly expressed 
in terms of the gyrocenter pull-back operator Tgy. The physical interpretation of the gyrocenter pull-back operation 
is therefore given in terms of the time integration of the guiding-center Vlasov equation over the fast gyromotion time 
scale. 

The gyrocenter pull-back also generates the standard decomposition of the perturbed particle Vlasov distribution 
in terms of its adiabatic and non-adiabatic parts as follows. Up to first order in e, we find 



T„ 



p ^ .aFo 1 dFo 



(80) 



where higher-order corrections have been neglected and the guiding-center Poisson bracket is now expressed in terms 
of the guiding-center coordinates (X, E, /x, 9) with the guiding-center energy E replacing the guiding-center parallel 
kinetic momentum pii. Next, we introduce the decomposition 



dFo ^ 
dE 



^1 = li^sc) + Gu (81) 



where the first term represents the adiabatic contribution to Fi and Gi represents its non-adiabatic contribution. 
Lastly, we use the identity ((26|l to obtain the relation 



(82) 



where /o = Tgc-Fb- Here, we note that the adiabatic contribution naturally separates into a particle part (involving 
(^) and a guiding-center part (involving Tgc {4>gc))- 



C. Gyrokinetic Vlasov Equation 

By combining the gyrocenter Vlasov equations (|74p and (I76p we obtain the nonlinear gyrokinetic Vlasov equation 

= ^ + {T. 5*^v},. (83) 



where the effective gyrocenter potential 



2n \] y^^'' n dt 

1 I I c^gy^g 



contains nonlinear ponderomotive corrections to the linear scalar potential. 
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We note that in standard gyrokinetic theory Q, the polarization-drift correction (involving dgy^gc /dt) is omitted 
and only the second term is retained in ([M)) . The gyrokinetic Vlasov equation ([M]) describes the time evolution of 
the gyroangle-independent gyrocenter Vlasov distribution F(X,p||,t; /l) in a 4 + 1 reduced phase space, where the 
gyrocenter Hamilton's equations (in a uniform magnetic field) are 



— ^ cb — . " — 

X = TJ|| b + — X V^'gy and P|| = -gb-V'I'gy, (85) 



and the gyrocenter magnetic moment JI is an invariant. 

A common approximation for the gyrokinetic Vlasov equation (I83p is obtained by writing it in truncated form as 
dgyF/dt = and then expressing, first, the gyrocenter Vlasov distribution function as F = Fq + Fi and, second, using 
the decomposition ((5T|) to obtain a gyrokinetic equation for the nonadiabatic part Gi. The electrostatic version of 
the Frieman-Chen gyrokinetic equation Q is thus obtained from the truncated equation 

dgyFi dgyFo 



dt dt 

= 9 ?5 ( ^ - ^ ) (0gc) - ^ X V(0ge) • Wo, (86) 



dFo / dgc d\ cb 

where the gyrocenter kinetic-energy coordinate E [see (|7ip ] is used instead of P|| and the background distribution Fq 
satisfies the guiding-center Vlasov equation dgcFo/dt = 0. Next, we introduce the decomposition ((8T|). where we write 



dt V dE J ^ dt dE 

to obtain the electrostatic Frieman-Chen gyrokinetic equation 

While this equation offers great simplicity for many practical applications, it also suffers from several deficiencies 
[2] (e.g., it lacks energy conservation when combined with the gyrokinetic Poisson equation) which limit its use in 
numerical simulations of electrostatic plasma turbulence. 



VI. GYROKINETIC POISSON EQUATION 



When the nonlinear gyrokinetic Vlasov equation (j83p is combined with the gyrokinetic version of the Poisson 
equation, we obtain a set of energy-conserving equations that can be used for numerical simulations of electrostatic 
plasma turbulence {3| . 

The gyrokinetic Poisson equation is expressed as a moment of the gyrocenter Vlasov distribution F through a 
sequence of phase-space transformations from particle to guiding-center to gyrocenter phase spaces: 



V 



2, 



47r 



Jd'zfS'^q J d'Z F6l, ^ q J d^Z F {Jg^SD 



where summation over particle species is implied, 5'^ = (5'^ (x — r) implies that only particles located at the field position 
X = r contribute to the scalar field (/)(r), and = T~j,^^^ = 5^{1^ + p^^ — r) is expressed in terms of the guiding-center 
gyroradius vector p^^. The last expression in (j88p involves the gyrocenter push- forward operation 



"•"gy- 



£{-5*1, 5}gc - ^{5*2, .g}gc - ^ l^"!, {Si, gjgcj 



to (5g(., we obtain 



where the generating scalar fields 5*1 and 5*2 are defined in (pOj) and ([55)1 . When the push-forward operator is applied 

/T-1,53 ) = (^3 ^ ^ 2 ^ 

\ gy "gc/ \"gc/ ^ 2^72 



(^|<i>gc, {I'gc, {51)} 



gc 



gc/ 



— e — 



1 ] 1 dgy^. 



'j'sc-H^Le^ ) + (89) 



o \ \ ^ dt 
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where 5^^ = S^^ — {S^^} = p^^ • V^^ + • • ■ denotes the gyroangle-dependent part of 5^^. We note that polarization 
effects enter into the gyrokinetic Poisson equation (|88|) through the term (J'^yS'^c) = (Jgy^^gc^'^))'- guiding-center 
polarization enters through the difference (T"^.^^^) — 5'^ while gyrocenter polarization enters through the difference 
(Tgy^Jg^) — ((5g(,}. We further note that (|89p may also be obtained from the functional derivative of the effective 
gyrocenter potential (|84p : 

^*gyM ^ {'^ly^D^ (90) 
so that the gyrokinetic Poisson equation ([55]) may be expressed in terms of the gyrokinetic variational principle [l|, Q 



5 



0. (91) 



The existence of a variational principle for the gyrokinetic Vlasov-Poisson and Vlasov-Maxwell equations allows us to 
compute exact conservation laws by Noether method Note that, according to the functional derivative (PO)) . the 

gyrocenter polarization effects (associated with the difference (Tgy^Jg^.) — {5^^)) require that quadratic nonlinearities 
in the electrostatic potential (p be retained in the effective gyrocenter potential ([84|) . 
Next, by introducing the gyrocenter gyroradius vector 

Pgy ^ Tg/ (X + pg,) - (X + pgj , (92) 

the push-forward expression ((89)) may be written as J^^^S^^ = (5'^(X + p^^ 4- Pgy — r). When expressed in terms of 
(S'l, 5*2, • ■ •); the gyrocenter gyroradius vector ([92]) becomes 

Pgy = - e {5*1, X + PgJ^^ - e2 15-2^ X + PgJ^^ 
+ - 1^1, X + PgJ^J^^ + •■• 

where the gyroangle-dependent part of the gyrocenter gyroradius vector ([93]) is 

Pgy = - e X}g^ + . . . = e A X V$gc + • • • , (94) 
while, up to second order, the gyroangle-indcpcndent part of the gyrocenter gyroradius vector ([M]) is 

Pgy = - ^ ({^1' Pgclge) + Y ({^1' {^1' X}gc} ) (95) 



gc/ 



gel 



where corrections of order fl ^ were kept. Expansion of the last term on the right side of the gyrokinetic Poisson 
equation in powers of Pgy yields the expression 

- = 1 J d'^pF - {<1 J d'pPgyF + ■••^ , (96) 

where we have ignored FLR effects in the first term on the right side (i.e., (^g^.) ~> S^), and the second term represents 
the polarization density. By keeping only terms of first order in e and lowest FLR order, where 

0gc ~ Pge • V(/) and $gc ~ - Pg(, • b X V0, 
the gyrorangle-independent gyrocenter gyroradius vector (|95p is expressed as 
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which involves the E x B velocity and the polarization-drift velocity. 

Because of the mass dependence appearing in (P7)) , ion polarization effects in the gyrokinetic Poisson equation ^ 

dominate over electron polarization. Lastly, in the standard nonlinear gyrokinetic formalism Q , the polarization-drift 
contribution appears at second order in e and is, therefore, omitted from the first-order gyrocenter gyroradius vector 
((97)) . The standard gyrokinetic Poisson equation (|96p thus yields the following relation between the electron (particle) 
density rte and the ion (gyrocenter) density Ui. 



I 



_B2/(47rmi 



Ii 



V0 
47r 



(98) 



where I_l = I — b b is the perpendicular unit matrix. In 
side of is zero] appears in the limit B"^ / {An mifii) -C 
condition 



, we note that the quasi- neutrality condition [i.e., the left 
In this limit, (j98p becomes the gyrokinetic quasi- neutrality 



erii 



Hi V I 



(99) 



which relates the electron (particle) density rig, the ion (gyrocenter) density n^, and the electrostatic potential 
(through the ion polarization density). It is important to note that the ion (gyrocenter) density Ui must be defined as 
the moment of the full ion gyrocenter Vlasov distribution Fi [i.e., it is a solution of the nonlinear gyrokinetic Vlasov 
equation ([55]) ] in order to conserve the global energy of the gyrokinetic Vlasov-Poisson equations [6, 7]. 



VII. SUMMARY 



The guiding-center recursive Vlasov (gcrV) method yields results that are identical to the gyrocenter Lie-transform 
(gyLt) method. The Lie-transform method, however, offers several computational advantages. First, instead of 
computing each gyrocenter variable individually (gcrV method), the derivation of gyrocenter variables by gyLt method 
involves a single function Sk at each order e'^' {k = 1, 2, ...) of the perturbation analysis. Moreover, we point out that, 
for most practical applications, the guiding-center and gyrocenter transformation can be kept separate since 

TgcTgyF + (eG?gy + G?g,) 9„F + 

Second, the explicit use of the guiding-center and gyrocenter pull-back and push-forward (Lie-transform) operators 
provides us with a simple interpretation of the recursive Vlasov method: the pull-back operator generates a fast- 
time-scale integration of the Vlasov dynamics while the push-forward operator represents the polarization dynamics 
in the gyrokinetic Poisson equation. Third, a self-consistent set of gyrokinetic Vlasov-Poisson equations is obtained 
by the Lie-transform method since it is derived from a variational principle (which guarantees the existence of exact 
conservation laws). Lastly, the Lie-transform method can easily be generalized to the fully electromagnetic case Q. 
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APPENDIX A: GUIDING-CENTER TRANSFORMATION 



In a nonuniform magnetic field (where es = p/ Lb denotes the dimensionless ratio of the characteristic gyroradius 
to the magnetic nonuniformity length scale), the guiding-center phase-space transformation is defined in terms of the 
first-order components 

= -Po = - {2^lB/m^^Y^ p, (Al) 
Gi" = (mc/e)^ (ai : Vb + b- V X b) - p|| Po • (b • Vb) , (A2) 

Gi = Po- U VlnB + b-Vb - p (a^ : + b-Vx bj , (A3) 
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= -Po-R + ^-VlnS + : Vb 



2fiB 



(A4) 



where we used the definitions for the gyroangie-dependent rotating unit vectors p = cos^ei — sin^?2 and -L = 
v±/|vx| = dp/d9, expressed in terms of the fixed unit vectors (ei,e2,h = ei X 62), so that the gyrogauge vector 
R = V± = Vei •?2 is gyroangle independent and the dyadic (traceless) tensors 



ai 



(ei?! -6262) sin 29 + (6162 +6261) cos 26 



and 32 = / ai d9 are gyroangle dependent. In (IAl|) - (IA4|) . p^^ = T^^x — X. — Pq + es Pi + ■ ■ ■ denotes the guiding- 
center gyroradius, where denotes the guiding-center push- forward operator (see below) and Pi is the first-order 
correction to the lowest-order gyroradius vector pg. 

In a nonuniform magnetic field, the guiding-center pull-back operator is expressed as 



Tgc = exp 



and the guiding-center push-forward Tg^.^ can be used to construct the guiding-center gyroradius vector 



Po + es 



1 „ 1 



(A5) 



(A6) 



where p^ • Vpg = - i (p^ • V In B) p^ - (p^ • Vb • pg) b - (pg • R) dp J 39. 

The Jacobian for the guiding-center transformation is expressed (up to first order in e^) as 

d 

{J G?) 



J. 



J - 



B 



eeV-iBpa) ~ eg B 



(dG^^^ dC; dG{ 



89 



B[l + esJ-b.Vxb) = Bl 



(A7) 



The guiding-center Poisson bracket, on the other hand, is expressed in terms of two arbitrary functions F and G as 

'dF dG dF dG\ n (OF dG OF dG" 



dw dt dt dw J B \ d9 dp, dp 89 



B*. Qpii Qpii 



cb 
9^ 



VF X VG , 



(A8) 



where B* = B + {p\\c/q) V X b and = b • B*. 

Lastly, the guiding-center pull-back operator (jASP and the guiding-center Jacobian (|A7P can be used to obtain the 
push-forward representation of the particle Vlasov-moment integral 



1x11 = j <fzJx5''{^-v)J^,F = j d''ZJ^,l-^^x5\y. + p^,-v)F, 



(A9) 



where x is an arbitrary function in particle phase space and we used the identity 



— -^gc Tg(. S. 



If we now expand the delta function (5^(X+Pgp— r) in powers of Pg^, and integrate by parts, we obtain the guiding-center 
push-forward representation 



IXII = llV^llgc - l|Pgc Vxllgc + 

which enables us to write particle fluid moments in terms of guiding-center fluid moments. 
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